1 



Perturbative Formulation of Pure Space-Like Axial Gauge QED 
with Infrared Divergences Regularized by Residual Gauge Fields 

Yuji Nakawaki and Gary McCartor* 

Division of Physics and Mathematics, Faculty of Engineering 
Setsunan University, Osaka 572-8508, Japan 
* Department of Physics, Southern Methodist University 
Dallas TX 15215, USA 

(Received ) 

We construct a new perturbative formulation of pure space- like axial gauge QED in which 
the inherent infrared divergences are regularized by residual gauge fields. For that purpose we 
perform our calculations in coordinates x M = (x + , x~ , x 1 ,x 2 ), where x + — x° sin 9 + x 3 cos 9 
and x~ = x° cos 9 — x 3 sin 9. A- = A cos 9 + A 3 sin 9 — n-A = is taken as the gauge fixing 
condition. That framework allows us to show that the residual gauge fields are not canonical 
variables in a pure space-like axial gauge formulation, although they are necessary ingredients 
of the system. To overcome that difficulty we first construct, in the region 0^9 < -j, a 
temporal gauge formulation in which x~ is taken as the evolution parameter and A+, Ai(i = 
f,2) and their canonical conjugate momenta are independent canonical variables so that 
we can quantize them by canonical methods. As a result we can obtain the commutation 
relations of the ^"-independent Nakanishi-Lautrup field B, which we extrapolate into the 
axial region ? < 9 < 5, where x + is taken as the evolution parameter. We obtain a 
formal solution of the gauge field equations with the relevant integration constants specified 
on the basis of the free gauge fields in the temporal formulation. Then we make use of 
the solution to specify the integration constants appearing when we solve the axial gauge 
constraint equations. It turns out that, because the residual gauge fields are multiplied by 
the inverse Laplace operator - d i +g u , which becomes singular in the axial region, they give 

rise to vanishing contributions to the equal a; + -time canonical commutation relations. For 
this reason, we can uniquely quantize the physical degrees of freedom and we do not obtain 
any extra terms in the equation for the Fermi field in spite of having introduced the residual 
gauge fields. We have also obtained the conserved Hamiltonian by integrating the canonical 
energy-momentum tensor T^v around a suitable contour. The Hamiltonian consists of the 
usual terms, with the residual gauge fields in A^ as regulators, plus the Hamiltonian for the 
residual gauge fields which is determined by integrating a density over the hyperplane x~ — 
constant. We show in detail that, in perturbation theory, infrared divergences resulting from 
the residual gauge fields cancel infrared divergences resulting from the physical parts of the 
gauge field. As a result we obtain the gauge field propagator prescribed by Mandelstam and 
Leibbrandt. By taking the limit 9— we can construct the light-cone formulation which is 
free from infrared difficulty. With that analysis complete, we perform a successful calculation 
of the one loop electron self energy, something not previously done in light-cone quantization 
and light-cone gauge. 



§1. Introduction 

Axial gauges, n^A^ = 0, specified by a constant vector have been used re- 
cently in spite of their lack of manifest Lorentz covariance. 1 ) One reason is that the 
Faddeev-Popov ghosts decouple in the axial gauge (AG) formulations. 2 ) For that 
reason we will not need to discuss the Faddeev-Popov ghosts in this paper and the 
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term ghost will always refer to residual gauge (RG) fields which are introduced as 
integration constants in solutions to constraint equations. The case n 2 = 0, the 
light-cone gauge, has been extensively used in the light-cone field theory in attempts 
to find nonperturbative solutions to QCD. 3 ) The same formulation is often used 
in phenomenological calculations where the low energy properties of the theory are 
parameterized in light-cone wave functions. Traditionally, constraints in AG for- 
mulations are solved to eliminate dependent fields in terms of physical fields. That 
elimination requires one to utilize inverse derivatives which introduce so-called spuri- 
ous singularities. It was first pointed out by Nakanishi 4 ) that there exists an intrinsic 
difficulty in the AG formulations, the solution to which requires an indefinite metric; 
even in QED. It was also noticed that in order to bring perturbative calculations done 
in the light-cone gauge into agreement with calculations done in covariant gauges, 
spurious singularities of the free gauge field propagator have to be regularized, not as 
principal values, but according to the Mandelstam-Leibbrandt (ML) prescription 5 ) 
in such a way that causality is preserved. Shortly afterwards, Bassetto et al. 6 ) found 
that the ML form of the propagator is obtained in the light-cone gauge canonical 
formalism in ordinary coordinates if one introduces a Lagrange multiplier field and 
its conjugate as RG degrees of freedom. 

It is now established 7 ^ that AG formulations are not ghost free, contrary to 
what was originally expected, and is still sometimes claimed. Nevertheless, the 
problem of finding the correct method by which to introduce RG fields into pure 
space-like AG formulations has not been completely solved. A first step toward 
solving this problem was made by McCartor and Robertson in the light-cone formu- 
lation of QED. 8 ) In previous papers, 9 ) we considered the solvable Schwinger model 
and learned ways to introduce and to quantize RG fields. In this paper we apply 
these methods to the more realistic problem of QED and construct an extended 
Hamiltonian formalism of QED in the ±-coordinates ), where 

x + = smQ + x 3 cos6>, x~ = x°cos6 — x 3 sin# and the space-like constant vec- 
tor n is specified to be = (n + , n~, n 1 , n 2 ) = (0,1,0,0) so that the gauge fixing 
condition is pure space-like. The same framework was used by others to analyze 
two-dimensional models. 10 ) We show that if we properly choose constituent fields, 
then we can construct the extended Hamiltonian formalism of QED in quite the 
same manner as the Schwinger model. 

In spite of the common use of light-cone gauge (with A + treated as a constrained 
field) and light-cone quantization, the one loop electron self energy has never been 
successfully calculated in that formulation. To our knowledge, that most funda- 
mental of loop calculations has been successfully calculated in three formulations of 
light-cone quantization: Morara and Soldati 7 ) calculated it in the gauge A + = 0; 
Langnau and Burkardt 11 ) calculated it in Feynman gauge but did not start the cal- 
culation on the light-cone; in 12 ) it was calculated in Feynman gauge starting from 
the light-cone and also calculated in light-cone gauge, but with a higher derivative 
regulator so that A + was a degree of freedom, not a constraint. But it has not been 
calculated in light-cone gauge with A + treated as a constraint. The problems that 
prevent such a calculation can be glimpsed in 11 ) and are discussed in detail in. 12 ) 
Basically, the infrared singularities are too strong to allow a successful calculation. 
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In the present paper we shall show that with the RG fields in place, the infrared 
singularities are softened in such a way as to allow a successful calculation. 

The paper is organized as follows. In §2, we construct the temporal gauge (TG) 
formulation and obtain the commutation relations of the Nakanishi-Lautrup field B, 
which are then extrapolated into the axial region. In §3 we consider the dipole ghost 
field and the constituent fields corresponding to those in the n-A = gauge Schwinger 
model. We then express A^ in terms of these constituent fields. In § 4 we obtain 
the translational generators by applying McCartor's method. 13 ** In § 5 perturbation 
theory is developed and we show that if we define the singularities resulting from 
inversion of a hyperbolic Laplace operator as principal values, the worst infrared 
divergences resulting from physical parts of the gauge fields are cancelled by infrared 
divergences from the RG parts. We then find that the remaining infrared divergences 
are regularized as the ML form of gauge field propagator. In that section we also 
perform a calculation of the one loop electron self energy. § 6 is devoted to concluding 
remarks. 

We use the following conventions: 
Greek indices [i,v, - ■ ■ will take the values +, — , 1, 2 and label the component of a 
given four-vector (or tensor) in the ± coordinates; 

Latin indices i,j,--- will take the values 1, 2 and label the 1, 2 component of a given 
four- vector (or tensor) in the ± coordinates; 

the Einstein convension of sum over repeated indices will be always used; 

x ± = (x ± ,x 1 ,x 2 ), x± = (x 1 ,^ 2 ), d 2 x± = dx 1 dx 2 , d 3 x = dx 1 dx 2 dx ± 
k± = (k±, hi, /C2), d 5 k± = dkidk2dk± 

g = cos 29, g ^ = g^ = sin 29, g ++ = — cos 29 

9-i = 9i- = 9+i = 9i+ = 0, gij = -5ij. 
7 + = 7 sin 9 + 7 s cos 9, 7" = 7 cos 9 — 7 s sin 9. 

§2. Temporal gauge formulation of n-A = gauge QED 

In this section we keep 9 in the temporal region 0^9 < ? and choose x~ as the 
evolution parameter. This enables us to show that RG fields are necessary ingredients 
in the canonical quantization of n-A = gauge QED. 

The n-A = gauge QED is defined by the Lagrangian 

L = --F^F^ - B(n-A) + V{i^D^ - m)V (2-1) 

where = d^ + ieA^ and B is the Lagrange multiplier field, that is, the Nakanishi- 
Lautrup field in noncovariant formulations. 14 ^ From ()2-l|) we derive the field equa- 
tions 

d„F» v = n u B + J v , J v = e&Y& (2-2) 



{i^D^ - m)& = 0, 



(2-3) 



4 



Y. Nakawaki and G. McCartor 



and the gauge fixing condition n-A = 0. The field equation of B, 

(n-d)B = d^B = 0, (2-4) 

is obtained by operating on (|2-2|) with d v . 

Canonical conjugate momenta are defined to be 



, SL 5L , SL 

7T+ = — — — = + , 7T~ = T7—T- = 0, 7T l = — — = F 



5d-A+ T ' <5d_A_ ' <59_A 



7rs = ^ = ' " = e = * 7 " (2 ' 5) 

and the Gaufi law constraint is described in terms of them as 

d + 7r + + di^ = B + J-. (2-6) 

We see from (|2-5|) and 1)2-6(1 that in the TG formulation we have three pairs of 
bosonic canonical variables, A + ,tt + , Aj,7r J (i = 1,2), which indicates indispensabil- 
ity of degrees of freedom other than the physical ones in constructing the canonical 
formulation. By using the canonical equal x~-time quantization conditions imposed 
on the independent canonical variables and the expression of B obtained from (|2-6|) 
we easily obtain the commutation relations of B: 

[B(x),A^y)] = id^(x + - y+), [B(x), vr + (y)] = [B(x),^(y)] = 0, 

[B(x),B(y)] = 0, [B{x),V{y)\ = e5^{x + - y + )V(y). (2-7) 

It is important to note here that, due to (|2-4|) . these commutation relations hold 
even when x~^y~ . 

The Hamiltonian, that is, the translational generator for the x~-direction, is 
given by 

P-= J d 3 x + TS = J d i x + {2-\^ + ) 2 + nj V " n+F +l ) 2 + n_(F +4 ) 2 + (F 12 ) 2 ) 
+ #(m - i 7 +9+ - Ydip + J^} (2-8) 

so that it is straightforward to develop an x~-time ordered perturbation theory. It 
was shown previously that the free gauge field is described as 15 ) 

Af_t = Tjj — -AB + d^A (2-9) 

V rp 

where is a free field satisfying 

□T M = 0, 0% = O, n-T = 0, (2-10) 
[T,(x),T„(y)] = i(- 9fiv + + _ n ^)D(x - y) (241) 

and 

V T 2 = + d 2 2 + n^d 2 , (2-12) 
A = ~(c- n_x~B - ^0d + B^j . (2-13) 
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Note in l|2-ll|) that D(x — y) is the commutation function of the free massless field 
and that the inverse derivative d~ is defined by 

1 f°° 

d- 1 f{x) = - dy-e(x--y-)f(x+,y-,x ± ) (2-14) 

^ J-oo 

which imposes, in effect, the principal value regularization. Note furthermore that 
C in (|2-13j) is the conjugate of B and satisfies 

d.C = 0, [C(x),C(y)]=0, [B(x),C(y)]=-iV T 2 S^(x+-y+). (2-15) 

It is useful to point out here that the physical degrees of freedom are carried 
by Tp, as is seen from (|2-lUj) . and the remaining degrees of freedom are carried by 
the pair of functions B and C. It should be noted in particular that solving the 
on-mass-shell condition 



= p 2 = n_p 2 + 2n + n_p + p_ — n_p 2 — p 2 , (p 2 = p 2 + p 2 ) (2T6) 



yields 



Ipl + n^pl - n + p + _ r— ^ n +yP+ + n -P± ~ P+ 

P- = , P =\/P+ + n-pl, p 



(2-17) 



from which follows 



1 n+P+ + yP+ + n-pl 
P- p 2 + ri-p? 



(2-18) 



Therefore we obtain the factors (p 2 + ri-p 2 ) 1 from the physical operators. It turns 
out that they are canceled by those resulting from — ^ , which is appled to B and C 

so that we obtain the following ML form of the x~-ordered gauge field propagator 
(0\T(A„(x)A u (y))\0) = J d'kD^k)^^ (2-19) 



where 



IT'U^v Tlj/hu 2 k^k u 



k z + is \ k- + 2£sgn(fc_|_) [k- + iesgn{k + )) z J 

We see from (|2-2()j) that RG fields are indispensable to regularize singularities asso- 
ciated with the gauge fixing in such a way that causality is preserved in complex fc_ 
coordinates. Because the factors (p 2 + n_p_j?) -1 , which becomes singular in the axial 
region due to n_ < 0, drop out completely from the gauge field propagator, we can 
formally extrapolate the gauge field propagator into the axial region. It is expected 
from this that we can develop a perturbation theory free from infrared divergences 
in the AG formulation. In next section we consider introducing RG fields so as to 
regularize the infrared divergences in the AG formulation. 
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§3. Constituent operators in the axial gauge formulation 

3.1. general solution of the gauge field equation \2- $Hj) 

We begin by obtaining the general solution of ()2-2|) in a way which is independent 
of the evolution parameter. Let a M (x) be a field satisfying 

□a„ = J„. (3-1) 

Then it holds that 



d„. „ . rip „ <9 (J 

t~2 

r 



^ - =AM< + ^ B = a M " v7^> ( V ± = ^ + 9 2 2 ) (3-2) 

V I V T> Vi 



which is verified as follows: If we apply D'Alembert's operator, which is described 
in our formulation as 

□ = n _d_ 2 + 2n + d+d- - n_d+ - V ± 2 , 

to both sides of (|3-2j) . then we obtain the same result, J M — ^^SjJj. Therefore the 

difference between the left hand side and the right hand side is a solution of the free 
massless D'Alembert's equation. The solutions which tend to zero at spatial infinity 
(|cc_|_|— >oo) vanish identically . 

In addition we can show that satisfies 

d^Of, = 0. (3-3) 

In fact operating on (J3-2|) with <9 M yields 

8^8 1 
- ^fdiAi = 3% - ^diJi 

.-. diipAi - d^A^ - Ji) = dii&'Fm - Ji) = -v ± 2 a^a M . 

The left hand side of the second line vanishes due to the spatial component of the 
gauge field equation and so does the right hand side. 
^From ()3-2|) is described as 



8 8 n 

A^ = a^- -^d m + -^=X - -Jt-B (3-4) 
v, Lv, 2 V T 



where we have denoted , 1 2 djAi as X. Throughout this paper the operator 



Vj 2 is understood as in the Fourier transforms of the relevant operators. 
The minus component of Q3-4JI has to vanish identically due to the gauge fixing 
condition. Imposing A- = gives rise to the following constraint 

°=X = -{a_ - J^aJ + ^B (3-5) 

-V, 2 V -L V T 
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the first term of which is rewritten, due to c^a^ = d + a + + d a_, as 
J^fl -_ 1 

v, 2 m ~ 

1 



a_ — -^-^diOi = — ^-^■{□a + — <9 + (<9 + a_ — <9_a + )} 



r{J_ - d+(d+a_ - 9_a+)}. (3-6) 



v ± 2 

Thus integrating (|3-5|) with respect to x~ yields X expressed as 

1 d- 1 

X = - — 2-{a+(9+a_ -9_a+) - J-} + A (3-7) 

V i 



where we have introduced integration constants involving C in the same manner as 
in (|2-13j) . Substituting (|3-7|) into X in (|3-4j) . we have expressed as 



71 

A^ = a^- - ^B + d„A. (3-8) 

It is straightforward to show that, because d-B = <9_C = 0, X and satisfy 
the following equations: 

8^8 8 

M 8-X = -n^B + -^diJ t -J-, (3-9) 

V i 



?7 77 i 

0%, = C - n_x~B + —^-d + B - (d-)- l J-. (3-10) 

We see from Q3-9J) that X is nothing but the dipole ghost field. 

3.2. Commutation relations of the constituent operators 

In the axial region, where ? < < ^ and x + is chosen as the evolution param- 
eter, we obtain from IJ2-1JI the canonical conjugate momenta 

*' = «^I = ' s^-*» + - (3 ' u> 

If we impose the gauge fixing condition, A^ = 0, on the equation for ir~ , then we 
obtain a constraint 

n~ = -<9_A+ (3-12) 
in addition to the Gaufi law constraint 

8-tt~ +dy = J+. (3-13) 

As a consequence, A + becomes a dependent variable so that there remain only 
Ai,ir l (i = 1,2) as independent bosonic canonical variables. Therefore, as quantiza- 
tion conditions we have only equal x + -time commutation and/or anticommutation 
relations on the independent canonical variables; the nonvanishing ones are 

[Mx),7r\y)] = iStjS^ix- - y-), {^(x),^(y) 7 +} = ^(aT - y~). (3-14) 
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The degrees of freedom of the system should not change when we move from 
the TG formulation to the AG one, because the field equations and the gauge fixing 
condition are the same. How can we resolve this paradox? The only answer is that 
B and C take the place of the canonical variables A+, tt + in the TG formulation but 
happen not to be canonical variables in the AG formulation. This reflects the fact 
that B and C are introduced as x~ independent static fields so that x + can not be 
the evolution parameter for them. Therefore, we cannot obtain their quantization 
conditions from the Dirac procedure. To supplement the insufficient quantiza- 
tion conditions, we assume that the commutation relations of B given in Q2-7J1 can 
be extrapolated into the axial region. This is a reasonable assumption, because we 
have introduced degrees of freedom for c-number residual gauge transformations and 
because B generates them. However it is not straightforward to make the extrapola- 
tion because the Laplace operator V T 2 becomes hyperbolic so that its inverse gives 
rise to singularities. Note that n_ = cos2# < in the axial region. We regularize 
these singularities by the principal value method. Consequently we obtain vanishing 
contributions at equal x + -time due to (V r 2 ) -1 . For example, we obtain 

[-Ls(x), -^=X{y)] x+=y+ = =^5^(x + - y + )\ x+=y+ 
1 jk ± -(x ± -y ± ) 



d * k+ ——^ e *"±-(*±-V±) = 0, (3-15) 



(2tt) 3 J + |fc ± | 2 + n_£^ 

: 1 T B(x),#(y)l + =y + = ^ 3 \x + -y+My)\ 3 



[ d s k + . u * u2 e ik ^-VMy) = (3-16) 



(27T) 

where we have made use of the fact that 



roc 



/-> wr^r ° (317) 

as a result of regularizing the singular integral by the principal value method. In 
this way B makes vanishing contributions to equal x + -time commutation relations. 

To express A + in terms of the independent canonical variables, we integrate 
(|3T2(I and (|3-13|) with respect to x~ . By integrating (|3-13|) we obtain 

vr~ = -d-A + = Tj-{J + ~ dy) + constant. (3-18) 

The integration constant is determined by comparing (|3-18|) with that given by A + 
in (|3-8j) and we get 

d+a- - d-a+ = —(J+ - din*), (3-19) 

nrr 

constant = ^d + B. (3-20) 
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Integrating (|3-18|) with respect to x~ and then comparing the result with A + in 
(|3-8|) . we obtain 

1 ■ 71 \ 

A + = -Q2( J+ - d '^) - ^2 B + 3 + A - ( 3 ' 21 ) 

We see from (|3-21[) that A + consists of the conventional physical operator plus 
RG operators. Furthermore, we see from (|3-7|) and (|3-19j) that <9+a_ — <9_a + plays 
fundamental roles in the axial gauge formulation. If we define U by 



d+a_ - 8-a+ = - drf) = ^-V ± 2 £, (3-22) 



then d + U is expressed as 

1 n -\/- V _L 2 1 

d + U = , d-djAj 1— 2 B + - J_, (3-23) 

-V ± 2 V t ^-V ± 2 

so that at x + = y + we obtain 

[E(x), d + S{y)} = i8® (aT - y~), [£(x), = = 0. (3-24) 

Here we have assumed the equal x + -time current commutation relations: 

[J+(x),J+(y)} = [J + ,J-(y)] = [«/-(*), J-(y)] = 0. (3-25) 

In addition, both U and d + U are gauge invariant so that due to Q2-7JI their four 
dimensional commutator with B vanishes: 

[E{x),B{y)\ = 0, [d + E{x),B{y)} = 0. (3-26) 

Note in particular that, when we calculate the commutation relations of d + S with 
X and we do not obtain any 5^ 3 \x + — y + ) contributions due to (|3-15|) and (|3-16|) . 
Consequently we obtain the following equal a; + -time commutation relations: 

[E{x),X(y)\ = -±-5^{x--y-), [d+Z(x), X(y)] = 0, (3-27) 

[Z(x),n i (y)] = 0, [d + Z{x)y(y)\ = ^=d-d^{x- - y~), (3-28) 

-V, 2 



[S(x),V(v)] = -e-@=LL6®(x- -y-)*(y), [Z(x),J-(y)] = 0, (3-29) 

-v 2 



[d+E(x), V(y)] = 6 ^— 6^ (a,- - y-) 7 +7_^(y), [d+S(x), J_(y)] = 0, (3-30) 
[S(x), J + (y)} = [d + E(x) } J + (y)} = 0. (3-31) 
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Now we can calculate the commutation relations of C, which is introduced as 
the integration constant in Q3-7|) . We rewrite ()3-7j) to express C in terms of operators 
whose commutation relations are known; we get 



2 



C = —E — J_ - X + -i =— n-X~B + f d+B . (3-32) 



Then using l3~T5|) . ([3~To]) and (l3~24l ~ (TOT) we obtain at x+ = y + 

[J^C(x),E(y)] = l^C(x),d + XJ(y)} = [^C(x), <%)] = 0, (3-33) 
r ' C(x),X(y)] = [-^C(x),J + (y)] = [-^C( 2; ),J_(y)]=0. (3-34) 



\/ rj-i \^ rp \/ rp 

If we use (j'2 • 7f) and (|3-26|) , then we obtain the following four dimensional commutation 
relation 



l —C{x) : B{y)] = [—^X(x),B(y)} =iS^(x+-y + ). (3-35) 



L v 2 

It follows from (|3~33l ~ (|3~35|l that 

[^C(x), ^Lc(y)]\ x+=y+ = 0. (3-36) 

It is clear that Af=Ai — -^-^Aj and 7r l ± =TT l — ^i-^ are gauge invariant and thus 

commute with B and C at equal x + -time. Therefore B and C do not give rise to any 
unwanted contributions to the equal x + -time canonical commutation relations. Now 
that we have Aj-, ir±, E, d + ZJ, B and C as fundamental operators, we can reconstruct 
our formulation in terms of them. To do this, we divide A^ into physical and RG 
parts as follows 

n. 



Ay = — ^B + d^A (3-37) 

where 



-V 2 



n ^ = ^(^-J + ), (3-38) 

T_ =0, Ti = A±+ jl |L I d + Z - jl J_ | , C!,-!9) 

0, I .. - 1 



7T 



+ ' , J+ . (3-40) 



rjT \ \l ^ rp" 
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It is easily seen that satisfies 

d»U„ = J v , {fftv = dfjTu — dvTft), (3-41) 
3% = — (3-42) 

and that Tj and tt 1 satisfy the equal x + -time canonical quantization conditions 

[T^Tjiy)] = [ir\x)^{y)\ = 0, [T\(x)rf{y)\ = »V (3) 0*~ - y~), 

[T i (x),V(y)} = [iT i (x),V(y)}=0. (3-43) 

§4. Translational generators in the axial gauge formulation 

In this section we start with the canonical energy-momentum tensor 

= -F vu d^A a + 9 ^F^F pa + M^d^ (4-1) 

and obtain the conserved translation generators. Note that in our formulation we 
have to resort to a nonstandard way of deriving the conserved generators. This 
is because A p in (|3-37|) has the term depending explicitly on x~ as well as x~- 
independent terms. As a result, the traditional formula for the Hamiltonian contains 
divergences. This reflects the fact that the integral d-T^~ dx~ does not vanish, 
although x~ is one of space coordinates. From the divergence equation, d v T^ = 0, 
we obtain 

ir/^ = 0, (4-2) 

where the integral is taken over a closed sur- 
face. From this we shall derive the conserved x ° 
generators. For the transverse directions we 

are justified in assuming that the integral x 1 

J^^diT^dx 1 vanishes. (Here, repeated in- 
dices do not imply sum over i.) Therefore, 
as the closed surface we use the one shown in 
Fig. 1, whose bounds T and L are taken to oo 
after the calculations are finished. A similar 
procedure was first used by one of the authors 
of the present paper. 13 ) It is straightforward 
to obtain 

= J d 2 x ± (^f L dx- [T+{x)]f=Z T + [T^( X )}l_Z L _^ (4-3) 

where 

[ T , + ( X )] X X + + Z-T = T+(x)\ x+=T - T/(x)U += _ T , 

[V0C-=-l = T »~( X )\*-=L ~ T-{x)\ x - = _ L . (4-4) 
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Let us illustrate the case of \x = + in detail. Because has the RG term, we 
divide T + + and into terms containing only physical operators, terms containing 
products of physical operators and RG operators and terms consisting solely of RG 
operators. Then, we rewrite these terms, which are to be integrated by parts, in the 
form of derivatives minus terms resulting from the integration by parts. As a result 
we obtain 

T + + = F+8+Ai - \{{F- + f + F+F +i + F~ t F-i - (F 12 ) 2 } + iW"y+d+9 

= \{{f- + ? + f^n - + (/i 2 ) 2 } + mi + d + * - dj+o+A 

+ J + (-^l B + 9 + A) - d_ (^/_ + T + + T + ^d + B + Ti—zdiB) 

I / 71 71 71 I \ 

+ di (/+(T+ + 8+A)) - - (^— 2 d + B^d + B + —d.B—^By (4-5) 

T + - = F_ + d + A + + F-fi+Ai + W^d+V 
= f'iUi + &-Y-D+9 + -L J" ^B - J ^B + B^d + C 

+ 5, (ri(T + + + d+Ti^-^B + ^B^^ + 2 B - ^B J-^+c) 

+ d + (f- + (T + + d + A) - om^B + ^ b^'b - ^*^M*) 

- a + + T 9+ + ^ B ^ B ) (4 " 6) 

When we substitute (|4-5p and (|4-6|) into ()4-3j) with /x = +, the derivative terms with 
respect to x % give vanishing contributions, whereas the last x~ independent term of 
(14-51) is cancelled out by the last linear term of x~ in (|4-6|) . We also see that B and 
C contained in the derivative terms with respect to x + in the fourth and fifth lines 
of (|4-6|) give vanishing contributions in the limit T— >oo because x + is one of spatial 
coordinates for B and C. As to B=f-^d+C in the second line, the contribution from 

the upper bound L cancels out that from the lower bound —L. However we keep it, 
because it is the Hamiltonian density for the RG fields. We see furthermore, that the 
derivative term with respect to x~ in the third line of (|4-5(l gives rise to a vanishing 
contribution, partly because its first term is cancelled out by the corresponding one 
in the derivative term with respect to x + in the fourth line of (|4-6j) and partly because 
the physical T + and T, vanish in the limit L^oo. 

With the redundant terms eliminated, we turn to consideration of the regular- 
ization term. The first term in the third line of (|4-5|) is the relevant one. We see 
that for this to work as the regularization term we have to get rid of —dif~^d+A, the 
remainder of the term F^d+Ai, in the second line of (|4-5|) . This is because if we 
combine this with J + d+A, then we get 

(J+ - dif^d+A = -a_/_+a+yl = d-(T + d-d + A - f-+d+A) (4-7) 
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so that all the derivative terms with respect to x~ in (j4-5[) are cancelled out by 
corresponding derivative terms with respect to x + in (|4-6|) . This results in losing the 
necessary regularizing term. We remark that Robertson and McCartor 8 ) overlooked 
this fact so that they did not obtain the regularization term in the light front for- 
mulation of QED. We have to discard the unnecessary term in such a way that the 
divergence equation is preserved. It turns out that if we supplement T + ~ with 



^4 



d. 



\j + d+d-A- ^J+d+A 



a. 



(4-8) 



and if we denote T + + with —dif^d + A subtracted as # + + , and T + with (|4-8j) added 
as 9_T, then we obtain the following divergence equation: 



d + e+ + a_0 H 

where 



djif+ifji) - di(f-+fl) - - di(Ji(T + + d+A)) (4- 



1 



{U-+) 2 + f+if + i-f-if-i + (fl2) 2 } 



+ \P(m - i~i~d- - i'fD^W + J + 
-l 



Ha 



V T 2 



B + d+A), 



1 



+ i—n_~$'y 7 (m — ij d_ — i~f l Di) , P 
n 



B-J- 



V T 2 



B + 0- [ ^J+8+d^A - ^- J+d + A 



(4-10) 



Ea. ([4-9[) indicates that it is enough to consider only # + + and . 

Now taking the limits L— >oo and T^oo in Q4-3JI gives us a conservation law 
which contains RG fields: 







d 3 x~ 



+ / d 3 x~ 



+ B 



1 



V T 2 



8:C 



(4-11) 



The derivative term in 1)4-8(1 gives a vanishing contribution to (|4-llj) and any other 
terms of are assumed to vanish in the limits x~^±oo in accordance with the 
traditions of the AG theories. Then we obtain that the first term is conserved by 
itself. We can add to it the constant Hamiltonian for the RG fields and thus 
we have 

(4-12) 



P + = Jd 3 x-0 + + (x) + Jd 3 x + B(x)^d + C(; 



We can similarly derive conservation laws for the — and transverse directions as 
follows: 



d 3 x~ [0_ 
d 3 x- [9, 



+ d a x 



2 V, 2 



-i x =oo 



B 



T 



+ / d 3 x' 



+ B 



1 



d,C 



— oo 



(4-13) 
(4-14) 
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where 



/V-i + »*7 + a-#> (4-15) 

\{(f- + f - f + J +l + Hf-i + (/12) 2 } + i*rd-V 

J-^B-O^J^B), (4-16) 

f+fiTj + i&i+diP, (4-17) 
f-+di + f-^Tj + #7-^ 
ft , 1_ „ 1 ^ „ , 1 



J -^2 B ~ J~diA + d + I —j+did-A-— J + diA) . (M8) 



From (|4-13jl and (|4-14|) we obtain the generators of translations for the — and trans- 
verse directions: 

P-= I d 3 x-9_ + (x) - - [ d 3 x + B(x)^B(x), (4-19) 
Pi = Jd 3 x-0 i + (x) + Jd 3 x + B(x)^- 1 d i C(x). (4-20) 

We end this section by showing that the Heisenberg equations for the constituent 
fields hold. For that purpose we use the expressions 

/-+ - ^-(ftvr - J ), f j = ' J+i = (_ n _^ » 7 i = (_ n _) » ( 4 " 21 ) 
to express the Hamiltonian in terms of Tj and 7r* as follows 
P+ = Jd 3 x-6 + + (x) + Jd 3 x + B(x)^d + C(x) 

= i jd 3 x-{{dl\dy - J+)) 2 + (-n_)" V - n+/_,) 2 + (-n.)!/^) 2 + (/ 12 ) 2 } 

+ [d 3 x-{&(m - i-f-d- - ifDi)& +J + (d + A - ^B)} + fd 3 x + B^C. (4-22) 

Note here that we can use the equal x + -time commutation relations to calculate the 
commutation relations with the operators contained in 9 + ^~(x), while B and C are 
regarded as independent of other constituent operators so that they commute with 
the others at all times. As a consequence we obtain 

f)- 7T 2 — 71 1 f 

[Ti(x),P+] = i^idrf - J + ) + i +p = i(d,T + + f +i ) = id+Ti, (4-23) 

[9(x),P+] = e(T+ - ^±*B + 8+AW - ^—(m - i-y~d- - i^D^ix) (4-24) 

where we have used the equality (7 + ) 2 = — ra_ as well as the fact that <? commutes 
with B and C contained in # + + at equal x + -time because they are multiplied by — 
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Equating 1)4- 24j) with d+\P and then multiplying both sides by j + provides us with 
the field equation of \P. 

In quite the same manner we obtain 

[T i (x),P r ] = id r T i , [V(x),P r ]=id r V. (r = -l,2) (4-25) 

Equations for B and C are calculated by the same rules and we get 

[ _L- B (x),P r ]=i^B, [A_ C (x) : P r ] = i^C, (r = +,l,2) (4-26) 

1 In 

[ v ^B(x),P^]=0, [^ C (x),P.] = -i^B. (4-27) 

It follows that 

[A l (x),P^]=id^A i . (4-28) 
§5. A perturbation formulation free from infrared divergences 

5.1. calculation of the propagator 

Now it is straightforward to develop an x + -time ordered perturbation theory by 
employing the following free and interaction Hamiltonians: 

H = \ J<£x-{{dl l dyf - nr V " n+f-if - n-U-if + (/12) 2 } 

+ jd 3 x~^(m - - 17%)$ + jd 3 x + B^C, (5-1) 

Hi = jd 3 x~{.PA^ - \j + ^J + ] (5-2) 

where we have denoted the free fields in the interaction representation with the same 
notation as those in the Heisenberg representation. Thus 

A^T^ + r^, r^ = --^B + d^A (5-3) 

V rp 

is the free gauge field. The physical part, T fl = — jpO-, is described in terms of 
free fields <9 + a_ — cLa + = t/-V 1 2 Z' and af = <n — as 



-V, 2 



a+ -?±a- = -±-—E, (5.4) 

g. q. g+ 

0, Ti = ai- -^-a_ = af H ■=J=—E. (5-5) 

d - /-V 2 d - 



If we consider the fact that the conjugate momentum of T, is given by 

vr* = f\ = d+af- + 9i d-E, (5-6) 

V, 2 
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then we see that E, d + U and af, d + af are two pairs of canonical variables. As a 
consequence we can express in terms of creation and annihilation operators as 

1 f rl^k - 2 

T " (i) = vm? J tsfE^h^*-)'-** + h.c. } (5.7) 



where 



— n + k + n + k- — k + 



k+ = \ k- ~ n - k2 \ , k = — , k+ = — , kj_ = \/ ki + k% 

V - 1 - n_ n_ V 

(5-8) 

and the operators a A (fc_) and a A (fc_) (A = 1,2) are normalized so as to satisfy the 
usual commutation relations, 

[a A (fc_), a y (q_)} = 0, [a A (fc_), a A ,(g_)] = <W<5 (3) (fc- - qr_). (5-9) 

Note that 

#>(*) = (0, 0. (5-11) 

are the polarization vectors satisfying 

^eW(fc) = 0, n^ A )(fc) = 0, (A = 1,2) (5-12) 

£ = + + - n 2 ^. (5-13) 

A=l 

We expand B and C in terms of zero-norm creation and annihilation operators 
as follows: 

1 1 f rl 3 k , 

~^B(x) = -j== J -jJ=9(k + ){B(k + )e-* k -* + £t(fc + )e^}U- =0 , (5-14) 



where 

[B(fc+), Ct(g + )] = [C(fc+), Bt (g+) ] = _ 5 (3) (fc+ _ (5 .i 6) 

and all other commutators are zero. We note here that limiting the integration 
region to be (0, oo) is indispensable to obtain the ML form of gauge field propagator. 
We define the vacuum state and physical space Vp, respectively by 

B(k+)\n) = C(k+)\n) = 0, (5-17) 



V P = { |phys) | £?(fc + )|phys> = }. 



(5-18) 
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Now we can calculate the x + -ordered gauge field propagator 
D^x -y) = (n\{6(x + - y+)A ll ( X )A u (y) + 9(y + - x^A^A^x)}^} 



1 



(27TY 



d^qDMe-^-ti. (5-19) 



It is straightforward to show that its physical part is given by 

p q z + ie \ q- q± J ql 

where q 2 = n_q 2 _ + 2n + q + q^ — n_q^ — q 2 ^_. We investigate in detail how the RG 
fields play roles as regulators. In the case that [i = i and v = j we obtain the RG 
contribution 

(nirCitfaO^fc,)) |Q> = -L^ ADf.(<?)e-^-a) (5-21) 

where 



/ r n - 

4( 9 ) = Mj j o & + p g _)^ 



/c i — Oi — ie /c+ + — 



2n_n + fc + / i z \ i 

" 5(g - } (n_^ + ^)2 \jfe + _ 9+ _ ie + fc + +g + -ie; J ' (5 " 22J 

Note that the explicit x~ dependence gives rise to the factor 5'(q-). Note also 
that there is no on- mass-shell condition among the RG field's momenta k + ,k\,k2 
so that there remains a fc+-integration. As a consequence there arise singularities 
resulting from the inverse of the hyperbolic Laplace operator. Nevertheless, when 
we regularize the singularities as the principal values, the integral on the first line of 
(|5-22|) turns out to be well-defined. In fact we can rewrite its integrand as a sum of 
simple poles and make use of the integral formulas 



OO 







dk+ ■ - — = -7rsgn(g+), (5-23) 

\k + — q + — ie k + + q + — ie 

P / dk + ( — i — ^ — ) = 0, (5-24) 



o 



k+ — a k+ + a 



where a = —M^=. As a result we obtain 



n_ / i i \ n_7rsgn(g + ) 



dk + T - - = 2 ta ^ ■ (5-25) 

n_k\ + q\ \ k + - q + - ie k + + q + - ie J n^q\ + q]_ 



On the other hand, the integral on the second line of (|5-22j) yields a linear 
divergence, as is seen by rewriting its integrand as a sum of simple and double poles: 

2n_i_n_A;+ 



+ 



(n-k\ + q\) 2 \k + - q + - ie k + + q + 



ie 



2n + n_(7 + / i i \ n + n^q\ — q\ 



(n-q\ + q^) 2 \k + - q + - ie k + + q + - ie J a {n_q\ + qj_) 2 \k+ - a 

i \ n+ ( i i \ , 

+ 77— 772 • ( 5 - 26 ) 



k + +a) n_q\ + qj_ \(k + - a) 2 (k + + a) 
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The integration of the first and second terms on the right hand side can be carried out 
with the help of (|5-23[) and (|5-24[) . However we cannot regularize a linear divergence 
resulting from a double pole by the principal value prescription. 17 - 1 We show below 
that this linear divergence is necessary to cancel a corresponding one that occurs in 
the physical part. For later convenience we rewrite the linearly diverging integration 
in the form 

00 / 1 1 \ f°° 2n(n^k 2 -q 2 ) 

Substituting these results and (|5-27f) into (|5-22j) yields 

g n-q iqj .n+qiqj f 00 2n_ (n_ k\ - q\ ) 

B%{q) = ^-5 (q-)n Sg n(u) ~ («-) / q ^ {nk 2 + + g 2 )2 , (5'28) 

where we have made use of the identity 

Thus, for the sum of (j5-20|) and (|5-28j) we obtain 

Aj(<?) = 2 , ■ ( ~ ftj - n 2 ^r- - m_g^ i 7rsgn(g + )(5 / (g_) 

. f 00 2n_(n_&l _ q 2\\ 
-n+qiqjd{q-) / ciA; + — 2 2 , 2 . (5-30) 

Now we can demonstrate that the linear divergence resulting from is canceled 

by the final term of (|5-3()|) when we restore Dij{x) by substituting (|5-3()|) into (|5T9|) . 
It should be noted here that q_ is conjugate to the spatial variable x~ , while k + is 
conjugate to the temporal variable x + . To make the infrared divergence cancellation 
mechanism work, both integration variables have to be either spatial or temporal. 
We show in Appendix A that if we change the integration variable from the spatial 

q_ to the temporal k + = n+q ~ — v^?- n -9 ±^ then the following integral formula holds: 

00 1^ 1,2 J2 



1 [<*> 2n+(n-k 2 + -ql) 



From (J5-3H) we see that the linearly diverging terms of the inverse Fourier transform 
of H5-3U|) are given by (we omit —in 2 qiqj for the moment) 

00 1/1 r°° 1 \ . _ 

00 g—iq—X y r°° \ 



-00 



q 2 . (q 2 + ie) n 2 q\ + q\ J q 2 _ 



^ ( ql q 2 +ie + q-{n 2 ql + ql){q 2 + ie) ) " (5 ' 32) 



-00 



Perturbative Formulation of Pure Space-Like Axial Gauge QED 



19 



We see from this that the last integrals diverge at most logarithmically, but logarith- 
mic divergences can be regularized by the principal value prescription so that there 
arise no divergences from (|5-32|) , This verifies that the following identity holds: 



~Y + i7rsgn(g + )5'(g_) - %_) 



Pf-9- + iirsgn(q + )5' (q^ 



1 



(q- + iesgn(q + )) 2 



(5-33) 



where Pf denotes Hadamard's finite part. Now substituting (|5-33|) into (|5-30|) yields 
the ML form of gauge field propagator: 



DM 



q 2 + ie 



-9ij 



n 2 qiqj 



(q- + zesgn(g + )) s 



(5-34) 



For other cases we omit detailed demonstrations because the calculations are 
similar. In the case that [i = + and v = i, we obtain the following RG contribution 



D g + M) = * 



dk A 





+ S'(q- 



n 1 



n_/q: + q L \k + — q + — ie k + + q + — ie 
ri-k + ( i i 



U-k\ + q\ \ k + — q + — ie k + + q^ 
2n-n + k\ ( i i 



l€ 



{n-k\ + q 2 ^) 2 \k + - q + - ie k + + q + - ie 
Because the integral on the first line gives rise to the imaginary part of 



we obtain 



(5-35) 



1 

<3_+jesgn(<j + ) ' 



D +i (q) 



n+qi 



n 2 q + qi 



q 2 + ie \q- + iesgn(q + ) ((/_ + iesgn(<? + )) 2 
For the case that jj, = u = + the RG contribution turns out to be 



(5-36) 



dk^ 



n_ 



_ k\ + qj_ V k + - q + - ie k + + q^ 



ie 



S(q-] 



2n + k + q 2 _ 

{n-k\ + qj_) 2 \k + - q + - ie ' k + + . 



ie 



(5-37) 



so that we obtain 



D 9 ++ (q) 



1 



2n + q + 5(q-)irsgn(q A 



q z + ie 

+n 2 q\5' (<?_)7rsgn(g + ) - q±S(q-) 



dq-- 



(5-38) 



Here it is noteworthy that the last term in (|5-38|) also cancels the linear divergence 
resulting from the contact term in ()5-2U(l . In fact it also yields a contact term, 
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as is seen from —5{q-) ^^_ ie = 5{q~) (l + ^q~^ J • Therefore, combining it with the 
corresponding one in (|5-2U|) and carrying out the inverse Fourier transform, we obtain 

i r «, U-) r - 4) = - rv^^i - — 



2vr7 V J -oo T- qlj o qi 2 

(5-39) 

Consequently, taking into account the fact that the Fourier transform of is 
— — + t — ^raOi we obtain 

i / 2n + g + n 2 g 2 



-C++ (?) = 2 l ■ -£++ + 



<? 2 + ie \ q- + iesgn(c/ + ) (<?_ + iesgn(g + )) 2 

This finishes the demonstration that, due to the RG fields, the linear divergences are 
eliminated even in the most singular component of x + -ordered gauge field propaga- 
tor. 

We end this subsection by investigating the light-front limit 9— ►? + 0. On the 
mass-shell, k 2 = 0, k + is given by k + = n+k ~ — V k -~ n - k ±_ go we g e j- 



( k 2 

lim k+ = I 2/t ( k - > °) (5.41) 
e-f +0 [00 (Jfe_ < 0) 

and contributions to in (|5-7j) from the integration region in which fc_ < vanish 
due to the Riemann-Lebesgue lemma. As a result, the integration region in 
is limited to (0,oo). It is worthwhile noting that in the limit, n_ = n 2 — >0, so we 
do not have the linear divergences except for the contact term in the most singular 
component of the gauge field propagator: 

£++(<?) = — ^—-\-i^5{q.){-i^gn{q + )) + ^5{q^ f°° dk+. (5-42) 
q- q 2 + ie q± qj_ q\_ J 

However, the linear divergence resulting from the contact term is also cancelled out 
due to an equality 

1 f°° „ 4 



dq- — - / dk + — = 0, (5-43) 
) q- J <7_|_ 

2 

which is verified by changing the integration variable from g_ to fc+ = J^. 15 ^ 

Zq— 

In the light-front formulation = -^j ^^ is a dependent field and has to 

be expressed in terms of the independent field tp + = -^7°7 + '!/>. We follow Morara 

and Soldati 7 ) to solve this problem and obtain the following free and interaction 
Hamiltonians 

H = jd 3 x-{±(f_ + ) 2 + i(/ 12 ) 2 + #7"5_ V} + jd 3 x + B^C, (5-44) 
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Hj = Jd 3 x-{J»A, + e^A^^-YA^ - \j + ^J + )- (5-45) 



The x + -ordered electron propagator is given by Chang and Yan 18 ) to be 



Sf = , n x 4 / d 4 p e 



(2tt 



-ip-x 



(p' + m) 1 7" 1 



p 2 — m 2 + ie 2 jr 1 



(5-46) 



We see here that there appears a noncovariant instant term, but its contributions to S 
matrices are cancelled by those resulting from the noncovariant fourpoint interaction 
term. 

5.2. calculation of the electron self energy 

In order to make comparison with previous calculations we shall regulate our 
calculation with the inclusion of Pauli-Villars fields. If we perform the calculation 
in the usual order for equal-time quantization, we should expect to need only one 
Pauli-Villars photon. However, it is common experience that more Pauli-Villars 
fields are needed when one performs the x + integral first. In Feynman gauge three 
Pauli-Villars photons 11 ) or one Pauli-Villars electron and one Pauli-Villars photon 
are needed. 12 ) In light-front gauge with higher derivatives so that A + is a degree of 
freedom, a much more complicated regulation is required including (in addition to 
the higher derivatives): three Pauli-Villars electrons, a Pauli-Villars photon and two 
carefully treated cutoffs. 12 ) In all these cases the Pauli-Villars electrons are included 
with flavor changing currents that break gauge invariance. The breaking of gauge 
invariance is removed by taking the masses of the Pauli-Villars electrons to infinity 
after the calculation is complete. It is not obvious that taking such a limit really 
restores gauge invariance but it can be shown to be true. 12 ) 

In the present calculation we keep one Pauli-Villars photon and investigate how 
many Pauli-Villars electrons are required to successfully regularize the one loop elec- 
tron self energy. We find that we need two Pauli-Villars electrons and so we use the 
Lagrangian 

1 1 1 

Y^(-\y +1 -F^ u F\ v -BA^+Y, -Mil^d^ - rmWi - ejrf^, (5-47) 

i=0 i 1 

where 

2 2 1 

^ = Z>' S> = ' ^ = E4- ( 5 - 48 ) 

i=0 i=0 i=0 

We shall have mo equal to the physical electron mass and vq = 1. For the other PV 
condition we shall require that {v{\ and {m^} be chosen such that the requirement 
from chiral symmetry — that the renormalized mass be zero if the bare mass is zero 
— be satisfied. That is, we shall require that 

<Mmo=o = 0. (5-49) 

Chiral transformations are dynamic in light-cone quantization and it is common 
that that requirement from chiral symmetry has to be imposed with an extra Pauli- 
Villars field. 18 ) As is usual, we shall regulate the infrared singularities with photon 
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masses and shall add these to regulate the propagator rather than add a mass term 
to regulate the Lagrangian. 

Using eqn's l|5-34|l . (|5-36j) . <|5-4U|) and 1)5-46(1 we find that the one loop electron 
self energy is given by 



(2tt)< 



/ A £(-i)^zv(g, /^KpItVI ( l P _ ft?" J 

i,3 1 



in; 



1 7 H 



+ ie 2 (p - q)+ 



(5-50) 

We shall take p± = 0. Due to the fact that they do not depend on masses, the 
contributions from the noncovariant term in the fermion propagator (the second 
term in Q5-46JI ). the contributions from the contact term in the gauge propagator 
(the last term in Q5-4UJI ) and the contributions from the second and third terms in 
(|5-45(l all go to zero in the sum over i and j. We can write the gauge propagator as 



"A 4 " g_+iesgn(g_|_) 

q 2 — fi 2 + ie 



(5-51) 



We then get 



(27T) 



We divide this expression into two pieces, U^\p) and U^ 2 \p), according to 



^ {2) (p) 



-ze 



(2vr) 4 

o2 



-IV 



(5-53) 



(2vr) 4 



J z — ' (p — q) — mj + is q z — \i ■ + ie 



We now perform the 7 algebra and take the inner product {plj S^\p)\p) to get 

,2 



ie 



„ „ 2£-(p + - o+) + 2^-(p_ - q) - Ami 



'J 



[(p - q) 2 - m i + * £ ] [<? 2 - m| + * £ ] 



Here it is important to note that when we perform the g+ integration, we obtain an 
extra contribution from the semicircle at infinity to the term propotional to (p + — q + ) 
as in the following 



dq A 



(p - q)- 



-ITT 



[(P ~~ l) 2 ~ m i + ie \ fe 2 - + i£ ] 

{sgn(g-) + sgn((p - <?)_)} ^+ m? 



4g- (p - c?)- 



P+ 



(p-q) 2 ± +mf~ie _ qj+vj- 
2(p-9)_ 2«j_ 



+ sgn(<?_) 
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and that it gives rise to a logarithmic divergence 



oo 



Aq_(p-q)_ 2p- J l-x 



so that a similar divergence resulting from 2(p-^q) — 1 * s cance Ued. As a consequence, 
changing the integration variable from q- to x = ^- yields 



< ™ <1, = (4/ rfV r & S^ ( - 1,i *(r 



V. 



mo(l — x) — 2m,j 

x) - mfx - fi?(l - x) — q\ 



I , r sr- . m (2x - 1) + mf - /xj 



By performing the x integration of the term in the second line we obtain 

mo(2x — 1) + — fi 2 



d q ± dx > i/^-iy— . g- - 

J Jo — rriQx{l — x) — mix — ^(1 — x) — q\ 

= - j d 2 q ± £ M-iy log = " / <*V £ „ log 



We see here that, due to the Pauli-Villars photon, this term is independent of Pauli- 
Villars electron masses so that if we impose the Pauli-Villars condition ^ Vi = 0, it 
is trivially zero. Consequently we obtain 



6m(1) = 7^3 / d ^ fdxY,n{-i) j — 

(2vr) 3 7 Jo Y mfr 



moil — x) — 2mi 



" x ) ~ m i x ~ - x) - «i ' 

(5-56) 

For terms involving only the physical electron mass we have 



i u 1 ^ 

The remaining contributions from are given by 

2 



2 r rl z 

PV (2vr) 3 y y± y y y m ° x(1 



mo(l — x) — 2mj 



x) - m?x - - x) - q\ 

(5-58) 

It is not difficult to show that this quantity goes to zero as the masses of the Pauli- 
Villars electrons go to infinity In this way we obtain the previous result that one 
Pauli-Villars electron and one Pauli-Villars photon are enough to regularize Sm^ . 

We now consider the contributions from U^ 2 \ Taking the matrix element, we 
find that 

" (2vr)4 mo] dq ^ Vt{ [q 2 - ^ + ie][q. + *£Sgn(g + )] 
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ie 2 f ,i s^s^ i ,» / 2(m -mi 



(27r) 4 7 "fr^ Up-^-mf + ieH^-^ + ie] 

2 — 1 J \ J 



2 m l- m 1 ~A 
mo ^ 



lip - q) 2 - m f + i £ ][q 2 - $ + ie][q- + i£sg n (<7+)] 



(5-59) 



The first term is trivially zero and the second term goes to zero as the masses of the 
Pauli-Villars electrons go to infinity. Thus we can discard them and concentrate on 
the third term, which is given by 

5m {2) = 5mf ] + 5m {2) (5-60) 

where 



-ie 2 f 2 [ dq- \- _ 
^+ f 

J —oo \ 



(p—q) 2 L +m 2 —ie q^+^—ie-i 
+ 2(p-q)- 2qZ \ 



bra. 



q_ f°° I 1 1 

* \q+-p++ 2(p x _ g) _ 
(2,_ *e 2 /\ 2 /• dg_ ^ 2^(-l)J^p+ 



l J ^+ 2(p- 9 )_ 2g_ -I 



IE 

X 



<Z 2 +e 2 ,_ 



/ I . (p-<?) 2 | +m 2 -ie <j 2 +M 2 -«£ I ' ^ ^ 

■ / -°° 2(p- g) _ ^-^—) 

By performing the g + integration Sm^ and <5m 2 2 ^ can be described, respectively, as 

2 2 

Sm? = f d 2 /" >_ m g +F (5 . 63) 

1 (2vr) 3 y y f-'g 2 +m 2 x + ^ 2 (l-x)-m 2 x(l-x) 

„2 



hJ 



( J m ° V ^ -(5-64) 



2p+Q-(P- -q-)~ q~{q± + m? - ie) - (p_ - + - ie) 



where 



q 2 ( gj + m 2 - 2ff + (p_ - g _)) 2 
(p_ + M;2) 



T = tan- ( , a e( tr/ t? "f f ?'7 g " )) J • (^6) 



Perturbative Formulation of Pure Space-Like Axial Gauge QED 



25 



Here we see that the term proportional to T is well-defined so that the factor g2 e +£ 2 

behaves as 7r<5(g_). As a result, it vanishes as e tends to zero. Therefore we can 
discard it. The remaining integral violates the chiral symmetry condition. That 
is why we need the second Pauli-Villars electron. We impose the last Pauli-Villars 
condition — Sm\ mo =o = 0. Actually, we see that 5m^ diverges as mo — > 0; so we 
set 



(2k) 



/"dvy>i(-i) j ' f dx- 

J i,3 J ° : 



x 



9 9 
i_ 

mo 



x* + e 2 q\ + mfx + /U 2 (l - x) 



+ 



7T 



m 



-p + log 



+ £ 2 Q- {q± + m i - «e) + (p- q)- (q± + - ie 



We can use this relation to write dm^ as 



0. (5-67) 



(5-68) 



where 



<5m 



(2) 
cl 



x 7725 — m. 
x' 2 + e 2 mo 

1 



qj_ + mfx + /x 2 (l - x) - infix (1 - x) qj_+ mfx + fi 2 (l - x) 



i,3 



x mj 



mr. 



x 2 + e 2 mo 



log 1 - 



mfix(l — x) 
mfx + /x 2 (l - x) 



(5-69) 



1,3 



ip+q-{p- -q-)- q-{q\ + n»? - *e) - (p- - + a*? - 



+ 2m (mo - mf)q-(p- - g_) log 



(p--g-)(gl + /x 2 ) 



{g_(gj + m 2 ) + (p_ - g-Kgj + ^ 2 )} 1 

2p+q-(p- -q-)- q-(q 2 ± + mf) - (p_ - g_)(gi + // 2 
,( 2 ) 



(5-70) 



Now the integrand of 5m c2 can be regarded as a continuous function of Conse- 
quently the factor 1 



2 , n behaves as <5(o_) so that we obtain 



lim <5m 



(2) 
c2 



(2vr) E 



mo 



log 1 



(5-71) 



26 



Y. Nakawaki and G. McCartor 



It follows from $EfM and (EZZH)) that 

lim flim5m (2)N ) = 0. (5-72) 

So the final answer for the electron self energy is regulated by one Pauli-Villars 
photon and is given by 

That is the same answer as the one obtained by Feynman methods. 12 )' 20 ) 

We end this section by making some remarks. If we regulate the spurious sin- 
gularity instead by 

* , (5-74) 
q- + iq + e 

which has been often used as the ML prescription, and if we take e to zero at the 
end of the calculation, then it improves properties of the q+ integration so that 
we need not introduce extra Pauli-Villars electrons. It is also possible to include a 
regulator in the spurious pole of the form +e ; in that case e can be taken to zero 
immediately after the q+ integration is performed and, if we keep e finite until after 
the Pauli-Villars electron masses are taken to infinity, we again obtain the Feynman 
answer. The main point is that the cancellation of the strongest singularity in the 
gauge propagator, shown in (|5-39() . allows a successful calculation of the electron self 
energy using standard techniques. 



§6. Concluding remarks 



In this paper we have constructed the quantization of QED in gauges fixed by 
specifying a constant, space-like vector, n, and the gauge condition, n^A^ = 0. We 
have then constructed the axial gauge formulation of QED in the zt-coordinates. 
Our framework has allowed us to consider axial gauges generally. The temporal and 
axial gauges in ordinary coordinates correspond, respectively, to 6 = and 6 = ^, 
while the light-cone formulation corresponds to 9 = ?. The most important aspect 
of this framework is that it has enabled us to use the temporal gauge formulation 
to obtain the algebra of the RG fields, which are not canonical variables in the pure 
space-like case; they are nevertheless, necessary ingredients in that case. 

We have obtained the commutation relations of B in the temporal region and 
extrapolated them into the axial region. We have also obtained the formal solution 
(13-81) of the gauge field equations and specified the relevant integration constants by 
comparing with the free gauge fields in the temporal formulation. We have made 
use of that solution to specify the integration constants that appear when we solve 
the constraint equations (|3-12() and 1)3- 13(1 in favor of the independent canonical vari- 
ables. We have specified the equal x + -time commutation relations of the constituent 
fields as being canonical and found eventually that the RG fields make vanishing con- 
tributions to the equal x + -time commutation relations because they are multiplied 
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by the operator We have furthermore obtained the conserved Hamiltonian to 

v T 

which Hamiltonian for the RG fields is added. In this way we have succeeded in 
constructing a perturbative formulation of pure space-like axial gauge QED in which 
the RG fields regularize infrared divergences inherent in the traditional quantizations 
of those gauges. The resulting gauge field propagator has the ML form. 

We have illustrated the effect of the RG fields by performing a calculation of the 
one loop self energy of the electron. In the usual light-cone gauge without higher 
derivatives so that A + is a constrained field, that calculation has not previously been 
performed successfully. The severe infrared divergences which result from solving the 
constraint equation for A + (without inclusion of the integration constants which we 
label B and C) have prevented a successful calculation. What we have found here is 
that including the integration constants (the RG fields) softens these severe infrared 
singularities and allows a successful calculation of the self energy using standard 
regulation techniques. The necessary regulation is slightly more complicated than in 
Feynman gauge but is considerably simpler than in the higher derivative regulated 
version of light-cone gauge where A + is a degree of freedom. Indeed we expect that 
the very complicated regulation procedures necessary in light-cone gauge without the 
RG fields 12 ) are somehow mimicking the effects of the RG fields, which should be 
included, to a sufficient degree that an effective renormalization can be performed. 
We should also remark that if another order of performing the integral is used (such 
as performing the q- integral first) then a smaller number of regulator fields is 
needed. We have studied the case of performing the q + integral first since that is 
the order that makes closest contact with nonperturbative light-cone calculations. 

While we have only considered the case of one loop (with the q + integral per- 
formed first), we expect that it may be possible to use the techniques of Paston 
and Franke 18 ) to show that the calculations with the relatively simple regulation are 
equivalent to Feynman methods to all orders. Since the quantization of QCD in 
Feynman gauge encounters difficulties which have not, so far, been solved, it may 
be that if the RG fields can be included in QCD, as they were here, the resulting 
formulation would have practical advantages over other formulations. 

Appendix A 

Verification of Eq. ()5-31(l 



We change the integration variable from g_ to k + = — ^ — q± . The quan- 
tity k + is a two- valued function of q- , and it takes its minimum value, yj— n_gj_=mo, 
at g_ = -^- m o- Therefore, when we change the integration variable from g_ 

to k + , we use q- = n+k+ ^+ m ° so that the region mo^k+ < oo corresponds 

to — co < q^^^—iriQ, whereas we use q- = n+k++ V k + m o go ^at the region 
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m,Q^k + < oo corresponds to ^r - m o=9- < °°- Hence we have 



d rco k + -n+Jk^-m , 
dk< = 



00 Q— J mo 



2 



(-n-)yk^ — m Q \ n + k + — \ k+ — m 



'0 



00 k + + n+\/ki. — m\ 



-n. 



+ / dk + ; ~ 1 , (A-i) 

Jm o (—nJ)Jk£ - m Q 2 \n + k + + Jk£-m 2 

where the first term diverges, but the second term is finite. Then, combining the 
first integral with the second one of (|5 31|) and making use of an equality 



(k + — k\ — m,Q)(n + k + + Jk~£ — m Q 2 ) 2 — 2n + ra_ \J k 2 — m Q 2 (n_ k\ — q\) 
= (k+ + n+y/kl - m 2 ){n+k + - ^ k 2 - m 2 ) 2 (A-2) 
yields 



°° dk k+-n+^h%-ml I _ n _ ^ + [°°dk 2n +( U - k + 



{-nJ)Jk* - m§ \n + k + -Jk*-mg) Jo (n_kl + q 2 ± ) 2 

00 k + -n + Jk 2 + -m 2 ( n + k + + Jk 2 -m 2 \ roc 2n + (n_/cl - q 2 , ) 
dk + ^^=^=- I r-h— — n I + / dk 4 



(-n-)yfh 



m o (—n-)\k_r — m, 



2_ m 2 \ n-k^ + ql I Jo (n_k 2 + q{) 



2 



k+ + n + Jk 2 + -ml I n+k + - Jk 2 -m 2 \ rm 2n + (n^k 2 - q 2 ) 
dk+ ; — I o \ + dk 



(-n.)^k 2 -m 2 \ J Jo {n-k 2 + q 2 ± ) 2 

(A-3) 

We see that the first integral in the last line is identical with the second one in (|A-1[) 
and thus is finite. Carrying out the integrations, we obtain 



f^^f'^ -"- = ^±, (A . 4) 

J™o (-n-)Jk£ - m Q 2 \n + k + + Jk£ - m Q 2 J n + q± 

; m '' dt+ ww) = _£EA. (A . 5) 

It follows that 

r dq^\ + r^+ 2 ? +( "/+~ 2 ?P = da + (us = 0. (a-6) 

J -00 q 2 - Jo {n_k 2 +q 2 Y 
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